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Abstract 

We study the electronic states of core multi-shell semiconductor nanowires, including the effect 
of strong magnetic fields. We show that the multi-shell overgrowth of a free-standing nanowire, 
together with the prismatic symmetry of the substrate, may induce quantum confinement of carriers 
in a set of quasi-lD quantum channels corresponding to the nanowire edges. Localization and 
inter-channel tunnel coupling are controlled by the curvature at the edges and the diameter of 
the underlying nanowire. We also show that a magnetic field may induce either Aharonov-Bohm 
oscillations of the energy levels in the axial configuration, or a dimensional transition of the quantum 
states from quasi-lD to Landau levels for fields normal to the axis. Explicit predictions are given 
for nanostructures based on GaAs, InAs, and InGaN with different symmetries. 
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Semiconductor nanowires (NWs) are emerging as new nano-devices with an enormous 
potential for applications. By now, several classes of materials can be grown in the form of 
single-crystal NWs, and key steps towards nano-electronic applications have been achieved: 
control of the position, monodisperse diameter, and uniformity of the NWai'^'^'^; n and 
p— doping, as well as growth of embedded barriers and superlattices along the growth 
direction;-!^!^ self-catalatio^'^'i^ or catalyst-free^i^iii growth which avoids deep traps forma- 
tion of extrinsic metal-assisted growthi^i^i^^; combination of different NW materials and 
substrates, including Si.- 

High quality free-standing NWs may also be used as substrates for epitaxial overgrowth 
of core shell structures. Also core multi-shell NWs have recently been demonstrated for 
GaN/In^-Gai^^N^^, GaAs/AlAs^^ and InP/InAsi^. Here, quasi-2D electronic systems can 
be formed, as in usual planar multi-layered heterostructures such as quantum wells, but 
wrapped around the NW. For example, in the InP/InAs structures^ electrons and holes are 
confined, due to the different band offset, in a few-mono layer shell of InAs overgrown on a InP 
core. Further developments in the growth and doping^^ may lead to an entirely new class of 
electronic devices and functionalities, based on high-mobility, bent two-dimensional electron 
gases with axial symmetry, confined in the low gap shells of the nanostructure. These may 
compete with carbon-nanotubes in the quest for the ultimate integration of nano-devices for 
solar cells^^ii^, opto-electronic applicationsi^i^^, sensing, or phase-based electronics. 

NWs typically grow with a prismatic shape, having the n-fold axial symmetry of the 
equivalent stable crystallographic directions. In fact, the stability of specific crystallo- 
graphic surfaces leads to faceting. Specifically, III-V semiconductor-based NWs (InP, InAs, 
GaAs)^'^'^ often grow as hexagonal crystals, exposing the {110} vertical facets normal to 
the [111] plane, but square (InP, InAs)2ii2^ or triangular (GaN)^i2^ shapes have also been 
demonstrated. Therefore, the electronic system confined in overgrown shells has the n-fold 
symmetry of the NW used as a substrate. For NWs with diameters in the range of few tens 
of nm, this results in a substantial deviation from the cylindrical shape which one would be 
tempted to use in modeling such systems. Indeed, one should note that a bent quantum well 
behaves like a weakly confining quantum wire in the region of the bent.^"^'^^ Therefore, the 
edges of a prismatic NW tend to be regions of preferred localization even with a constant 
well width, this effect being concurrent or in competition with possible width fluctuations 
in the region of the bent. As we show below, the nanostructure may behave as a set of 
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n tunnel-coupled quantum channels, depending on structural parameters at the nanoscale, 
such as the NW diameter and symmetry or the curvature at the edges. 

In this Letter we investigate the electronic states of core multi-shell NWs. Taking ex- 
plicitly into account their prismatic symmetry, we show that quantum confinement occurs 
in a set of quasi-lD channels located at the edges of the NW. An axial magnetic field af- 
fects the overall phase of these states and induce Aharonov-Bohm (AB) oscillations of the 
multiplets arising from the prismatic symmetry. A transverse magnetic field, instead, may 
deplete these channels and induce Landau level formation along the NW facets. Results are 
presented for several devices based on different classes of materials and symmetries. 

As we mentioned, a bent quantum well behaves as a ID electron waveguide in the area of 
the bent, since there is a larger region to accommodate the wave-function.— i21i2^ This local- 
ization effect does not require, and would add to, any fluctuation in the quantum well width 
(which is likely to happen in real samples^^). Solving the appropriate 3D Schrodinger equa- 
tion of the nanostructure would be numerically inconvenient in the present large systems, 
even within the envelope function approximation. In fact, we are focusing on electronic 
states which are conflned in a layer which is much thinner than the overall large diameter 
of the NW. Under these assumptions, the core of the NW acts as substrate for the epitax- 
ial overgrowth, determining the symmetry of the shells. This allows us to model the bent 
surface as a strictly 2D curved structure (see Fig. [1]). The conflned carriers move under the 
action of an effective potential which can be derived by a limiting procedure, starting from 
the full 3D effective-mass Schrodinger equation^"^ and taking the thickness of the bent 
layer down to zero. Such curvature-induced effective attractive potential 

Vhent = —- 1 2' V-*-) 

where m* is the effective mass of the carriers in the shell. The potential depends on the 
local radius of curvature of the bent. Note that due to the limiting procedure, electrons 
sit in a 5-like shell; therefore, the only parameters appearing in the above equation and in 
the system model presented in Fig. [1] are the shell curvature and radius and the effective 
mass of the carrier corresponding to the shell material. 

TEM images^ show that each edge of the NW actually consists of a region extending 
over a few nm with a flnite curvature. Accordingly, we model the cross-section of a prismatic 
NW starting from a polygon (see Fig. [1]) with a circumscribed circle of radius R, and we 
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FIG. 1: Left: the cross-section of a prismatic NW is modeled from a regular polygon (dashed line). 
The radius of the circumscribed circle is R. Each edge is rounded for a length / and with a constant 
radius of curvature Tc = al (see text). A transverse magnetic field B± forms an angle 6 with a 
line normal to a facet. The azimuthal position on the rounded polygon (solid line) is given by the 
X coordinate, at a distance d{x) from the center. 7(x) is the angle between the tangent in x and 
the line perpendicular to d{x), and s(x) is the distance between x and the direction of the field 
through the center of the polygon. Right: charge density of the ground state of a carrier in a core 
multi-shell NW (gray scale), showing preferred localization along the edges (darker regions). The 
state shown here is for the GaAs-based nanostructure of Fig. EJ^b). 

round off each corner for a length / at each side of the edge. Each rounded edge has a 
constant curvature l/vc, where = al, with a = y/S for the hexagon, a = 1 for the square 
and a = for the triangle. A spatial coordinate x is defined which runs on the rounded 

polygon (solid line in Fig. [1]), so that Vc = rc{x). The prism is therefore modeled as a n-fold 
succession of planar facets (where Vbent = since Tc oo) and rounded edges (with a 
constant attractive potential due to the constant value of Vc). Here we suppose an abrupt 
change of the curvature between the edges and the prism sides, that is not likely to happen 
on the atomic scale. Nevertheless, a progressive change between the two curvature would 
only influence the shape of the well defined by Eq. ([1]). The results presented in the following 
would remain qualitatively unchanged. 

Since ([T]) does not contain the wave-vector ky along the axial direction, the energy levels 
split in a set of parabolic ID subbands En + ffky/2m*, where the ky = subband edges En 
are simply determined by a ID effective mass equation in the step-like potential Ves{x) = 
Vhentirdx)). We assume that the curvature of the shell is different form zero only at the 
rounded edges. The resulting effective potential has a series of square wells (whose depth is 
given by Eq. ([T])) located at the edges, and separated by regions with Ves = corresponding 
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FIG. 2: (a) Lowest energy levels in the GaAs layer of an hexagonal core multi-shell NW with 
R = 40 nm at zero magnetic field, as a function of I, i.e. of the curvature at the edges, (b) 
Confining potential V^g (x) (bottom panel) and charge densities of the five lowest energy levels En 
corresponding to / = 2 nm (vertical dashed line of part (a)). Note that the first and second excited 
levels are doubly degenerate. 

to the prism flat sides. In Fig. [2]^a) we show the lowest energy levels En of an electron 
confined in a GaAs3^ shell overgrown on the surface of an hexagonal NW with i? = 40 
nm, as a function of /. The energy levels are easily understood as the result of the six- 
fold attractive potential felt by the carriers, shown in the bottom panel of Fig. [2](b) for the 
realistic valued I = 2 nm. The potential clearly resembles a multiple quantum well ~10 meV 
deep, with periodic boundary conditions. As shown by the charge density, the ground state 
of the system consists of a set of tunnel- coupled channels located along the edges of the NW, 
with a confinement energy in the few meV range. The probability densities corresponding 
to the four lowest energy levels show preferred localization at the edges, where the potential 
is attractive, but with a substantial probability density along the facets of the NW. 

For larger / the curvature diminishes, and the quantum wells become wider and shallower, 
until, for / = R/2 (for an hexagonal symmetry), a cylindrical system is obtained, where Vbent 
just amounts to an addictive constant. At the opposite limit of small /, the wells become 
thinner and deeper. The width and depth changes do not compensate, and the confined 
states depend non trivially on I (or, equivalently, on the curvature). At the envelope function 
level used here, an hexagonal NW belongs to the D^h space symmetry group. Accordingly, 
energy levels group in six-fold bunches, separated by gaps that increase for decreasing /, 
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which controls tunnehng between the edges. The ground state is not degenerate, while the 
first and second excited levels are doubly degenerate. The third excited state is again not 
degenerate and has six nodes in the zero-potential regions. The topmost graph of Fig. MJo) 
shows the lowest state of the second bunch, which has maxima located at the centers of 
the facets of the NW and has nodes in each well. In a cylindrical system the 3rd and 4th 
levels would be degenerated^ but the presence of the wells breaks the rotational invariance 
around the NW axis and lifts the degeneracy of those levels. Clearly, this behavior is 
exclusively determined by the periodic boundary conditions and the geometric symmetry 
of the system. The issues described above for a hexagonal wire can be easily extended to 
square or triangular NWs. 

We now consider the effects of a magnetic field. In this case the eigenvalue problem 
is more complicated, but a proper choice of the gauged provides an analytical expression 
of the 2D effective-mass Schrodinger equation, where the dynamics on the surface is still 
decoupled from the transverse one {i.e. normal to the curved surface). We first consider 
the magnetic levels of a semiconductor shell wrapped around a prismatic NW in an axial 
magnetic field of intensity B^\. As in the zero field case, the longitudinal and azimuthal 
directions are uncoupled, and the axial energy dispersion is parabolic. The ky = energy 
levels. En, are the eigenvalues of the effective ID Hamiltonian^^ for a charge q 

2" 

jd(x) c,os'y(x)B» 1 

2 



2m* 



-h'^dl + iqhB\\d{x) cos^{x)dx + ( -qd{x) cos'y{x)B\ 



Veffix), (2) 



which describes the quantum dynamics in the azimuthal direction. Here, d{x) is the distance 
of X from the center of the polygon, 7(x) is the angle defined in Fig. [H The familiar cylin- 
drical case d{x) cos7(x) = R is recovered when / spans one half of the polygon side.— >^>^ 

From Eq. ([2]) we have calculated the energy levels En as a function of the field intensity 
for a number of materials and geometries taken from the literature. In Fig. [3](a) we show 



the magnetic levels in the GaAs s 
R = 40 nm and / = 2 nm (see Ref. 



lell overgrown on the surface of an hexagonal NW with 
isl). Note that with an axial field only discrete rotations 
around the axis are symmetry operations,— and the magnetic symmetry group is lowered to 
Cgy. Accordingly, all energy levels are non degenerate. Due to the 6-fold symmetry of the 
system, however, levels form a series of braids containing six levels, each one separated from 
the higher one by an energy gap that depends on /, as shown in Fig. [2l In Fig. [3](a) dashed 
lines indicate the values of the field generating a flux $ multiple of the flux quantum $o- 
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FIG. 3: Magnetic levels En vs axial field intensity, Sy, at fcy = for a set of core multi-shell 
NWs. The values of corresponding to multiples of $o are indicated by vertical dashed lines, (a) 
Hexagonal GaAs shell with i? = 40 nm and / = 2 nm;— (b) Hexagonal InAs shell with R = 32 nm 



and / = 5 nm; (c) Square InAs shell (see Ref. |2ll | for the InP core) with R = 70 nm and I = 5 
nm; (d) Triangular Ino.25Gao.75N shell with R = 18.5 nm and / = 4 nm.— 

The energy levels show clear oscillations with a period $/$0; a periodicity that is typical of 
AB oscillations, as shown in cylindrical systems.— 

In Fig. [3]^b) we show the magnetic levels in the InAs shell^ with hexagonal symmetry^ 
with R = 32 nm and I = 5 nm, showing AB oscillations similar to the previous GaAs case. 
However, here the separation between the braids is smaller since l/R is larger (see the SEM 



images in Ref. 39|]), which enhances the tunneling between different edges of the structure. 
To highlight the effect of the different symmetry, (Fig. [3](c)) we calculated the magnetic 
levels for a hypothetical system, where a low gap layer of InAs is overgrown on the surface 
of a InP core NW with square symmetry.^^ The parameters are R = 70 nm and / = 5 nm 
(see SEM images in Ref. 2l|]). AB oscillations are still present, but, according to the square 
symmetry of the system, each braid is formed by four levels. 

Finally, in Fig. [3]^d), we show the magnetic levels of an Ino.25Gao.75N shell^^ with a 
triangular symmetry, with R = 18.5 nm and / = 4 nm.— As expected, AB oscillations 
appear and each braid contains three levels. We underline the different energy scale and 
field ranges in the four cases, which depend on the different values of R (influencing the 
energies En of the azimuthal modes and the magnetic flux and the value of m* that 
rescales the energies. 
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We next analyze the effect of a homogeneous magnetic field which is perpendicular to the 
NW axis with intensity B±. As in the previous case, it is still possible to write an effective 
ID Hamiltonian^ for the azimuthal dynamics, but now it also depends on the wave-vector 
ky along the axis, 

Hb± = ^ [-h'dl + {hky - qB^s{x)f] + V,ff{x). (3) 

Here, s{x) is the distance defined in Fig. [TJ Given the polygonal cross-section, the results 
will be also dependent on the orientation 6 of the field with respect to the facets. 

Since the component of perpendicular to the surface varies from facet to facet (de- 
pending on 9), there is no single magnetic length scale characterizing the system, contrary 
to the usual planar 2D electron gas case. As a consequence, the carrier quantum dynamics 
exhibits two different regimes, that can be analyzed starting from Eq. ([3]). The Hamiltonian 
contains a linear term in B±ky, and a term proportional to Bj_. The former dominates at low 
field or large longitudinal wave-vector and induces orbital Zeeman splitting.— It represents 
the classical Lorentz force, which induces carrier localization in a quasi- ID channel on one 
side of the NW, where the vertical component of the field tends to vanish.^^ For sufficiently 
high fields, the dominating quadratic term i?^ induces carriers localization in two strips 
located at opposite sides of the NW, where the field is perpendicular to the surface, these 
states being reminiscent of Landau levels.— These localization mechanisms compete with, 
or add to the localization along the edges of the NW, depending on the field orientation. 

By solving Eq. ([3]), we have computed the density of states (DOS) shown in Fig. HI for the 
GaAs (top) and InAs (bottom) shell with hexagonal symmetry. Similar results can be easily 
inferred for the square or triangular systems discussed above. For each system, the DOS 
has been calculated with two limiting orientations of the magnetic field: ^ = and 6 = n/G 
(see Fig. [T]). At each magnetic field, the magnetic energy levels form a set of subbands with 
a minimum (darker regions in Fig. H]) which, in general, corresponds to a state with ky ^ 0, 
due to the linear term in B±ky. The lowest subbands (below the levels labeled as Landau 
states) correspond to states which tend to be localized on one side of the NW with respect to 
the field direction, that is where the orthogonal component of B± on the surface of the NW 
vanishes. For ^ = (right panels) the localization area coincides with an edge of the NW. 
Therefore, edge-induced localization is enhanced by the field in one of the side edges of the 
NW, while the other edges are charge depleted. For 9 = n/G (left panels), on the other hand, 
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FIG. 4: DOS (gray scale) vs energy, as a function of the intensity of the perpendicular magnetic 
field B±. Top: hexagonal GaAS with R = 40 nm and I = 2 nm; Bottom: hexagonal InAs with 
R = 32 nm and I = 5 nm. Left panels: field oriented at ^ = ir/d. Right panels: field oriented at 
6 = 0. (see Fig. [T]for the definition of 6). 

the area of field-induced localization coincides with one of the NW facets, and all edges are 
charge depleted. Note the differences between the GaAs and InAs results, which are due to 
the different tunnel-coupling between the edge wells. At higher energies, peaks labeled as 
Landau states shift linearly with B^. These states are confined by the field on the top and 
the bottom of the NW with respect to the field direction. Again, if these regions correspond 
to edges of the NW (that is for 6 = 7r/6), localization is enhanced for these two edges, while 
the other four edges are charge depleted. On the other hand, if 6* = carriers tend to be 
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confined in the center of tlie two facets where the field is perpendicular, and all edges of the 
NW are charge depleted. Note the different slopes of the energy levels with field intensity 
for different orientation 9 of the field. In fact, only the orthogonal component cos(^)5x is 
effective, and the Landau-level energy is correspondingly smaller. For higher energies, the 
succession of peaks is repeated. A discussion of these states in a related cylindrical systems 
is given in Ref. 35|. 

In conclusion, we have shown that in core multi-shell NWs quantum confinement occurs 
in a set of quasi-lD channels located at the edges of the NW, and we have evaluated such 
confinement for typical parameters. An axial magnetic field induces AB oscillations of the 
multiplets arising from the prismatic symmetry. A transverse magnetic field may enhance 
edge localization or deplete these channels and induce Landau level formation along the 
facets or edges of the NW, depending on the field direction with respect to the NW facets. 

This work has been partially supported by project FIRB-RBIN04EY74 and Cineca Cal- 
colo Parallelo 2008. 
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